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like transcription factor 

Jerome S. Menet, 1 ' 2 Stefan Pescatore, 
and Michael Rosbash 2 

Howard Hughes Medical Institute, National Center for 
Behavioral Genomics, Department of Biology, Brandeis 
University, Waltham, Massachusetts 02454, USA 

The mammalian circadian clock relies on the master 
genes CLOCK and BMAL1 to drive rhythmic gene ex- 
pression and regulate biological functions under circadian 
control. Here we show that rhythmic CLOCKrBMALl 
DNA binding promotes rhythmic chromatin opening. 
Mechanisms include CLOCKrBMALl binding to nucleo- 
somes and rhythmic chromatin modification; e.g., in- 
corporation of the histone variant H2A.Z. This rhythmic 
chromatin remodeling mediates the rhythmic binding of 
other transcription factors adjacent to CLOCKrBMALl, 
suggesting that the activity of these other transcription 
factors contributes to the genome-wide CLOCKrBMALl 
heterogeneous transcriptional output. These data there- 
fore indicate that the clock regulation of transcription 
relies on the rhythmic regulation of chromatin accessi- 
bility and suggest that the concept of pioneer function 
extends to acute gene regulation. 
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Circadian clocks drive the rhythmic expression of a large 
fraction of the transcriptome in many eukaryotic tis- 
sues and organisms to regulate biochemical, physiologi- 
cal, and behavioral functions. Circadian gene expression 
is generated by a set of core clock genes, which interact 
in transcriptional feedback loops. In mammals, they 
include the two master heterodimeric transcription fac- 
tors CLOCK and BMAL1. This heterodimer rhythmically 
activates the expression of their transcriptional repres- 
sors, Period (Perl and Per2) and Cryptochrome (Cryl and 
Cry2) (for review, see Mohawk et al. 2012). Although 
these proteins and other clock components are well 
characterized, recent evidence suggest that the mecha- 
nisms by which they control genome-wide rhythmic gene 
expression are not well understood (Rey et al. 2011; 
Menet et al. 2012). 

We recently characterized CLOCK:BMALl target genes 
and analyzed rhythmic transcription using Nascent-seq in 

[Keywords: circadian rhythms; regulation of transcription; nucleosome 
positioning; nucleosome occupancy; MNase-seq; chromatin modifica- 
tions,- H2A.Z] 

1 Present address: Department of Biology, Center for Biological Clocks 
Research, Texas A&M University, College Station, TX 77843, USA. 
Corresponding authors 
E-mail rosbash@brandeis.edu 
E-mail menet@bio.tamu.edu 

Article is online at http://www.genesdev.org/cgi/doi/10.1101/gad.228536.113. 



mouse livers (Menet et al. 2012). The data revealed a 
surprising disconnect between the phase of CLOCK: 
BMAL1 DNA binding and the phase of target gene 
transcription, including the transcription of key core clock 
genes. CLOCK:BMALl binding occurs at the same phase 
of the cycle for all target genes, whereas the peaks of cy- 
cling transcription are heterogeneous, with little or no 
relationship to the singular phase of CLOCK:BMALl 
binding (Menet et al. 2012). This indicates that other 
transcription factors are involved in core clock gene 
transcription, for which there is experimental support 
(Ukai-Tadenuma et al. 201 1 ). 

There is evidence that CLOCK:BMALl function ex- 
tends beyond transcriptional activation. For example, 
CLOCK is reported to have histone acetyl transferase 
activity (Doi et al. 2006), and genome-wide rhythmic 
modifications of CLOCK target gene chromatin have 
been recently described (Feng and Lazar 2012; Koike 
et al. 2012; Le Martelot et al. 2012; Vollmers et al. 2012; 
Aguilar-Arnal and Sassone-Corsi 2013). Importantly, ec- 
topic expression of Drosophila CLOCK generates ectopic 
circadian clocks (Zhao et al. 2003; Kilman and Allada 
2009), indicating that CLOCK can function to establish a 
circadian program analogous to the developmental pro- 
grams generated by key factors like Pax6 (Osumi et al. 
2008). Because some transcription factors such as the 
glucocorticoid receptor (Truss et al. 1995; Nagaich et al. 
2004; Voss et al. 2011) establish their programs by remod- 
eling chromatin and promoting nucleosome removal at 
lineage-specific target genes (Magnani et al. 2011; Zaret 
and Carroll 2011), we hypothesized that CLOCK:BMALl 
plays a similar role at its target gene chromatin. 

In the present study, we report that the rhythmic 
binding of CLOCK:BMALl on DNA promotes the rhyth- 
mic removal of nucleosomes at its binding sites. Relevant 
mechanisms include CLOCK:BMALl binding to nucleo- 
somes as well as rhythmic chromatin modifications such 
as incorporation of the histone variant H2A.Z. We also 
show that this rhythmic chromatin opening at CLOCK: 
BMAL1 DNA-binding sites is associated with rhythmic 
binding of another transcription factor (HNF6). This 
suggests that the activity of other transcription factors 
contributes to the heterogeneous transcriptional output 
of CLOCK:BMALl target genes and that the activity of 
these other factors relies on the rhythmic regulation of 
chromatin accessibility of CLOCK:BMALl. 

Results and Discussion 

To test the hypothesis that CLOCK:BMALl remodels 
chromatin, we first performed a genome-wide nucleo- 
some analysis using MNase-seq (digestion of chromatin 
with micrococcal nuclease and high-throughput sequenc- 
ing of mononucleosomes) in mouse livers at six time 
points across the light:dark cycle. Four wild-type animals 
were used for each time point, and a minimum of 84 
million nucleosomes was sequenced per time point. 
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CLOCK:BMALl is a pioneer-like transcription factor 



Analysis of CLOCK:BMALl DNA-binding sites (3217 
peaks from Koike et al. 2012) showed that the nucleo- 
some signal is rhythmic at these sites for all four in- 
dividual rhythms and lower when CLOCK:BMALl binds 
to DNA during the light phase (e.g., Zeitgeber time 02 
[ZT02] to ZT10) (Fig. 1A ; Supplemental Fig. 1 ; Rey et al. 
201 1; Koike et al. 2012). In contrast, the identical MNase-seq 
assay in Bmall~^~ livers showed a more limited decrease 
in nucleosome signal and no rhythmicity, indicating that 
CLOCK:BMALl directly contributes to nucleosome re- 
moval (Fig. 1 A; Supplemental Fig. 2). Consistent with this 



Top 400 CLK:BMAL1 peaks 



B 




ZT02 

ZT06 

ZT10 

ZT14 

ZT18 

ZT22 
— Bmall -/- 



-0.6 -0.4 -0.2 0 0.2 0.4 0.6 
Distance from CLK:BMAL1 peak center (kb) 



Top 400 REV-ERBa peaks 



o0 0 oooo 



O Wild-type ZT6/10 

• Wild-type ZT18/22 
O Bmall -/- ZT6/10 

• Bmall -/- ZT18/22 



CLK:BMAL1 DNA binding strength 



Gene Body (top25%, n=406 peaks) 




ZT02 

— ZT06 
ZT10 

ZT14 
ZT18 

ZT22 
Bmall -/- 



-1 -0.5 0 0.5 1 

Distance from Rev-erbA peak center (kb) 



Intergenic regions 
(top25%, n=180 peaks) 




-0.8 -0.4 0 0.4 0.8 

Distance from CLK:BMAL1 peak center (kb) 



TSS (top25%, n=146 peaks) 




-0.8 -0.4 0 0.4 0. 

Distance from CLK:BMAL1 peak center (kb) 



■0.8 -0.4 0 0.4 0.8 

Distance from CLK:BMAL1 peak center (kb) 



Figure 1. CLOCK:BMALl promotes the rhythmic removal of 
nucleosomes at its DNA-binding sites. (A) Average nucleosome 
signal at the top 400 CLOCK:BMALl DNA-binding sites (±0.6 kb) 
in mouse livers during the light phase (ZT2, ZT6, and ZT10; green) 
and dark phase (ZT14, ZT18, and ZT22; red/orange) of wild-type 
mice and in Bmall" ^ mice (average signal for six time points; 
black). (B) Effect of CLOCK:BMALl DNA-binding strength on the 
average nucleosome signal at CLOCK:BMALl DNA-binding sites 
(±75 bp), at times of high binding (ZT6 and ZT10; open circles) or 
low binding (ZT18 and ZT22; closed circles) to DNA in wild-type 
(black) and Bmall" h (blue) mice. The 3217 CLOCK:BMALl peaks 
were equally distributed in 10 bins based on the strength of 
CLOCK:BMALl DNA binding. (C) Average nucleosome signal at 
the top 400 Rev-erbA DNA-binding sites (±1 kb) in mouse livers. 
[D-F] Average nucleosome signal at the top 25% of CLOCK:BMALl 
DNA-binding sites (±0.8 kb) located in gene bodies [D), intergenic 
regions [E) t and TSSs (F) in mouse livers. 



conclusion, the magnitude of the loss of nucleosome 
signal correlates with the strength of CLOCK:BMALl 
DNA binding (Fig. IB). 

We next examined nucleosome signal at REV-ERBa 
DNA-binding sites. This protein is a circadian transcrip- 
tion factor, and its gene is a direct CLOCK:BMALl target. 
Like CLOCK:BMALl, REV-ERBa rhythmically binds to 
DNA, with higher binding at the end of the light phase 
(Feng et al. 2011; Bugge et al. 2012; Cho et al. 2012). 
In contrast to CLOCK:BMALl -binding sites, however, 
REV-ERBa sites had indistinguishable nucleosome sig- 
nals at all time points. Moreover, the signals were not 
affected in Bmall~^~ livers despite a blunting of Rev-erba 
transcription and expression in this genetic background 
(Fig. 1C; Supplemental Fig. 3; Kornmann et al. 2007; 
Menet et al. 2012). These results show that not all rhyth- 
mic transcription factors promote nucleosome removal 
and suggest that CLOCK:BMALl transcriptional regula- 
tion is special. 

Interestingly, CLOCK:BMALl -dependent nucleosome 
removal is more pronounced within gene bodies and 
intergenic regions than at transcription start sites (TSSs) 
(Fig. 1D-F; Supplemental Figs. 4, 5). Consistent with this 
observation are the nucleosome signals from Bmall~^~ 
liver chromatin, which are strongly dependent on the 
genomic location of the CLOCK:BMALl peaks: TSSs are 
depleted of nucleosomes similar to wild-type chromatin, 
whereas nucleosome depletion in gene bodies and inter- 
genic regions is much less pronounced in Bmall~^~ mice. 
The low TSS nucleosome signal even in the absence of 
CLOCK:BMALl binding (e.g., during the night at ZT18 
and ZT22 as well as in Bmall~^~ mice) is consistent with 
the literature,- i.e., TSSs are generally more nucleosome- 
depleted because of intrinsic DNA sequence bias and/or 
the presence of the transcription machinery (Hughes 
et al. 2012; Iyer 2012; Thurman et al. 2012). In contrast, 
nucleosome removal at enhancers (within genes or at 
intergenic regions) is more dependent on cis-regulatory 
mechanisms such as sequence-specific transcription 
factors (Calo and Wysocka 2013). The data show that 
CLOCK:BMALl promotes nucleosome removal and opens 
chromatin more potently at enhancers than at TSSs. 

If a major role of CLOCK:BMALl is to open chromatin, 
how is this achieved? Nucleosomes present a physical 
barrier to most transcription factors and inhibit consen- 
sus sequence recognition (Magnani et al. 2011; Zaret and 
Carroll 2011; Dunham et al. 2012; Thurman et al. 2012). 
However, some developmental transcription factors, 
called pioneer factors, can bind to DNA within nucleo- 
somes and promote histone repositioning and/or removal 
(Magnani et al. 2011; Zaret and Carroll 2011). To address 
whether CLOCK:BMALl has similar properties, we per- 
formed a CLOCK chromatin immunoprecipitation (ChIP) 
on mouse livers DNA digested by MNase at two different 
time points, ZT06 and ZT22 (Fig. 2; Supplemental Fig. 6). 
These time points are the times of maximal CLOCK: 
BMAL1 DNA binding and when CLOCK:BMALl initi- 
ates DNA binding, respectively (Rey et al. 2011; Koike 
et al. 2012); importantly, nucleosome signals at CLOCK: 
BMAL1 sites are still maximal at ZT22 (Fig. 1). 

We found that CLOCK could immunoprecipitate 
mononucleosomes (Illumina libraries were size-selected 
to ensure an insert length corresponding to only one 
nucleosome) (see the Materials and Methods for more 
details), and the signal strikingly resembles those usually 
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ZT22 



CLK ChlP-Seq 
Input (ZT22) 




ZT06 



■ CLK ChlP-Seq 

■ Input (ZT06) 



0 0.5 1 -1 -0.5 0 

Distance from CLK:BMAL1 peak center (kb) 



MNase-treated chromatin 




Sonicated chromatin 



ZT06 
ZT22 




-1 



-0.5 



Distance from CLK:BMAL1 peak center (kb) 

Figure 2. CLOCK binds to DNA wrapped around nucleosomes. [A) 
CLOCK ChlP-seq signal on mononucleosome (i.e., mouse liver 
chromatin digested by MNase) at ZT22 (light blue; left) and ZT06 
(dark blue; right) for the top 400 CLOCK:BMALl DNA-binding sites. 
The signal corresponds to the average of three independent ChlP-seq 
experiments. Input MNase-seq signal at the same binding sites is 
displayed for both time points (ZT22 [dark orange] and ZT06 [green]). 
(B) CLOCK ChlP-seq over input signal ratio on MNase-treated 
chromatin at ZT22 (light blue) and ZT06 (dark blue) at the top 400 
CLOCK:BMALl DNA-binding sites. (C) CLOCK ChlP-seq over 
input signal ratio on sonicated chromatin at ZT22 (orange) and 
ZT06 (green) at the top 400 CLOCK:BMALl DNA-binding sites. 



The H2A.Z signal is strongly rhythmic at CLOCK: 
BMAL1 DNA-binding sites: high from ZT02 to ZT14 and 
then decreasing during the night to reach a trough at 
ZT22 (Fig. 3A). Importantly, incorporation of H2A.Z is 
severely compromised in Bmall~^~ mice, and levels do 
not exceed the trough levels observed in wild-type mice 
at ZT22 (Fig. 3 A; Supplemental Fig. 8). Furthermore, 
there is a striking correlation between H2A.Z signal and 
CLOCK:BMALl -mediated decrease in nucleosome signal 
(Fig. 3B). In Bmall~/~ mice, a higher H2A.Z signal at TSSs 
is associated with a stronger decrease in nucleosome 
signal compared with intergenic regions or gene bodies 
(Figs. IF, 3B). Conversely, a higher amplitude of H2A.Z 
signal occurs in wild- type mice within intergenic regions 
and gene bodies, which is associated with a stronger effect 



7 Intergenic regions (n=717) 

6 




-1.5 -1 -0.5 0 0.5 1 1.5 
Distance from CLK:BMAL1 peak center (kb) 



seen for histone modifications,- e.g., two peaks equidis- 
tant from the DNA-binding site and separated by —400 
base pairs (bp) (Fig. 2A). Importantly, CLOCK was asso- 
ciated with mononucleosomes at both ZT22 and ZT06, 
with only slightly higher signal/input ratio at ZT06 (Fig. 
2B ; Supplemental Fig. 7). To validate this result, we also 
performed a control CLOCK ChIP on the same mouse 
liver nuclei but sonicated rather than MNase-treated. As 
previously described (Rey et al. 2011; Koike et al. 2012), 
there is much stronger CLOCK binding to sonicated 
DNA at ZT06 than at ZT22 (Fig. 2C). The data therefore 
indicate that CLOCK:BMALl first binds to its target sites 
within nucleosomes (e.g., ZT22) and then promotes 
nucleosome removal to effect subsequent binding to 
naked DNA (e.g., ZT06). 

Nucleosome dynamics include histone post-transla- 
tional modifications and histone variants. Modifications 
often affect positively charged lysines and weaken histo- 
ne :DNA interactions (Cosgrove et al. 2004; Henikoff 
2008), and the H2A variant H2A.Z similarly weakens 
histone:DNA interactions and aids nucleosome reposi- 
tioning by chromatin remodelers (Henikoff 2008; Jin et al. 
2009; Ku et al. 2012; Li et al. 2012; Hu et al. 2013). To 
determine whether CLOCK:BMALl -mediated nucleosome 
remodeling involves H2A.Z, we performed a H2A.Z ChIP 
from mouse liver chromatin across a light:dark cycle. 
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Figure 3. Rhythmic CLOCK binding on DNA is associated with 
rhythmic H2A.Z signal at CLOCK:BMALl DNA-binding sites. [A) 
H2A.Z ChlP-seq over input signal ratio on MNase-treated chroma- 
tin in wild-type mice during the light phase (green) and dark phase 
(orange/red) and in Bmall~l~ mice (average of six time points; 
black). Signal ratio is displayed at CLOCK:BMALl peaks located 
within gene bodies, intergenic regions, or TSSs. (B) Average H2A.Z 
ChlP-seq/input signal (top) and percentage of nucleosome signal 
decrease [bottom) at TSSs, intergenic regions, and gene bodies. [Left) 
Values in Bmall~ h mice. [Right) CLOCK:BMALl -specific contri- 
bution (e.g., values above those observed in Bmall ^ mice). 
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of CLOCK:BMALl on nucleosome removal (Figs. 1D,E, 
3B). This is similar to what has been shown recently for 
FoxA2 during differentiation of embryonic stem cells into 
endoderm/hepatic progenitors (Li et al. 2012); namely, 
CLOCK:BMALl binding to nucleosomes promotes use of 
the histone variant H2A.Z and nucleosome removal. 

Nucleosomes at CLOCK:BMALl DNA-binding sites 
have been recently shown to be rhythmically acetylated 
and methylated (H3K9ac, H3H27ac, H3K4mel, and 
H3K4me3) (Koike et al. 2012). Interestingly, the phases 
of H3K27ac and H3K4me3 rhythms are influenced by the 
genomic location of CLOCK:BMALl sites (Supplemental 
Fig. 9) and better match the phase of CLOCK:BMALl 
DNA binding at gene bodies and intergenic regions (from 
ZT0.4 to ZT12.5) than the later phase of DNA binding at 
TSSs (H3K27ac phase: ZT16.6; H3K4me3 phase: ZT14.9). 
As the gene bodies and intergenic sites are also the most 
potent sites of nucleosome removal (Fig. 1D-F), the data 
suggest that these rhythmic histone modifications part- 
ner with H2A.Z incorporation to help promote nucleo- 
some removal, as proposed in other systems (Ku et al. 
2012; Hu et al. 2013). 

A simple model to explain the disconnect between the 
uniform phase of nucleosome removal at CLOCK: 
BMAL1 sites and the heterogeneous phases of transcrip- 
tion is that other transcription factors bind to the open 
chromatin at these sites. To address this possibility, we 
first analyzed 31 previously published ChlP-seq mouse 



liver data sets to identify transcription factors that bind 
close to CLOCK:BMALl sites (see the Materials and 
Methods). There is a significant coassociation of most 
transcription factors,- i.e., transcription factor DNA-bind- 
ing sites cluster in DNase I-hypersensitive sites (Fig. 4 A; 
as described in Dunham et al. 2012; Thurman et al. 2012). 
This analysis also confirmed that nuclear receptors 
associate better with the circadian repressors Cryl and 
Cry 2 than with other clock components (Lamia et al. 
2011; Koike et al. 2012). To extend this analysis, we also 
assayed the percentage of base-pair overlap between these 
31 mouse liver ChlP-seq data sets. This analysis revealed 
that many liver transcription factors bind adjacent to 
CLOCK:BMALl DNA-binding sites. Some factors even 
overlap with >70% of CLOCK:BMALl sites,- e.g., HNF4A, 
Bcl6, STAT5, CEBPA, or HNF6 (Fig. 4B). 

To compare transcription factor-binding dynamics ad- 
jacent to CLOCK:BMALl sites with nucleosome occu- 
pancy rhythms, we performed HNF6 ChIP assays at two 
opposite time points. Consistent with the model, HNF6 
rhythmically associates with DNA-binding sites located 
close to a CLOCK:BMALl DNA-binding site (Fig. 4C). 
The amplitude of binding was surprisingly high (greater 
than fourfold for many sites), suggesting that HNF6 af- 
fects CLOCK:BMALl transcriptional output. Importantly 
rhythmic binding does not occur at control sites,- e.g., 
HNF6 sites without nearby CLOCK:BMALl -binding sites 
(Fig. 4D). Assuming that the activity of other transcription 
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Figure 4. CLOCK:BMALl -mediated rhythmic nucleosome removal promotes the rhythmic binding of transcription factors to DNA. (A) 
Coassociation between transcription factors in mouse livers. Thirty-one publicly available mouse liver ChlP-seq data sets were analyzed by pairs 
using the Genome Structure Correction statistic as previously described (Dunham et al. 2012). Black rectangles denote core clock genes and 
nuclear receptors (see the text for more details). (B) Percentage of overlap between 31 publicly available mouse liver ChlP-seq data sets. Black 
rectangles denote transcription factors that exhibit an overlap superior to 40% with core clock genes (see the text for more details). (C) Rhythmic 
nucleosome signal at a CLOCK:BMALl DNA-binding site located near HNF6 DNA-binding sites. Nucleosome signal is displayed for wild-type 
mice at time of high (average ZT6 and ZT10; green) or low (average ZT18 and ZT22 ; red) CLOCK:BMALl DNA binding. The HNF6 ChlP-seq 
signal from Faure et al. (2012) is shown in gray. Genomic locations of CLOCK:BMALl (blue) and HNF6 (black) consensus sequences are also 
displayed. (D) HNF6 ChlP-seq signal in mouse livers at ZT08 (white) and ZT20 (black) at several specific DNA-binding sites [n = 4 mice per time 
point). Values represent the average ± SEM. (**) P < 0.05; (***) P < 0.01. [E] Model illustrating a new mechanism by which CLOCK:BMALl 
regulates the expression of its target genes. 
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factors is similarly potentiated, their nature (e.g., acti- 
vator or repressor) as well as their level of expression 
must contribute to the heterogeneity of CLOCK:BMALl 
target gene transcription quantitatively and qualitatively 
(Fig. 4E). 

Taken together, our data indicate that CLOCK:BMALl 
functions like pioneer transcription factors and regulates 
the DNA accessibility of other transcription factors. 
Contrary to the permanent chromatin opening associated 
with lineage commitment, however, CLOCK:BMALl- 
mediated chromatin opening is dynamic and occurs every 
day. This may reflect differences between the down- 
stream physiologically relevant transcription factors like 
HNF6 and the downstream developmentally relevant 
factors that follow pioneer proteins. 

This close parallel between CLOCK:BMALl and pio- 
neer transcription factors can explain why ectopic ex- 
pression of CLOCK in Drosophila leads to the develop- 
ment of functional ectopic clocks (Zhao et al. 2003; 
Kilman and Allada 2009): The chromatin of key core 
clock genes is opened like the chromatin of key de- 
velopmental target genes. It is notable that a large num- 
ber of additional CLOCK:BMALl target genes are proba- 
bly not core clock genes but output genes. Many of them 
also undergo comparable nucleosome changes but are 
only marginally expressed. CLOCK:BMALl may there- 
fore promote a rhythmically permissive state to accom- 
modate circumstances not yet encountered, also reflect- 
ing the absence of one or more transcription factors under 
these conditions. Last, the suggested pioneer mechanism 
for CLOCK:BMALl provides a different perspective on 
the core circadian clock: Circadian feedback may directly 
control the temporal regulation of key target gene chro- 
matin and only indirectly impact the transcriptional 
activation of many key core clock genes. 



Materials and methods 

Animals 

Adult wild-type and Bmall mice (Bunger et al. 2000) entrained on a 12 h- 
light:12 h-dark schedule were used. All experiments were performed in 
accordance with the National Institutes of Health Guide for the Care and 
Use of Laboratory Animals and approved by the Brandeis Institutional 
Animal Care and Use Committee (IACUQ protocols no. 0809-03 and no. 
12013). 



Generation of MNase-seq libraries 

Formaldehyde-cross-linked mouse liver nuclei, collected as previously 
described (Menet et al. 2012), were incubated with 1000 U/mL micrococ- 
cal nuclease for 15 min at 37°C before the reaction was stopped by 
addition of EDTA. These conditions resulted in ~70%-80% mononu- 
cleosomes and ~20%-30% dinucleosomes. Sequencing libraries were 
generated using 50 ng of DNA purified from the MNase-digested chro- 
matin (Illumina TruSeq DNA sample prep kit) and size-selected to ensure 
an insert size of a mononucleosome. Four MNase-seq libraries per time 
point (except for Bmall' 1 ' ZT22, for which n = 3) were sequenced. Each 
library corresponds to one mouse. 

ChIP 

MNase-digested chromatin was immunoprecipitated using H2A.Z (Ac- 
tive Motif, 39113) or CLOCK antibody (Santa Cruz Biotechnology, sc- 
6927X). The number of mice used was as follows: for H2A.Z, four wild- 
type and two Bmall' 1 ' mice per time point (six time points), and for 
CLOCK, three wild-type mice per time point (ZT06 and ZT22) 



Sonicated chromatin (Diagenode Bioruptor Plus) was immunoprecipi- 
tated using CLOCK or HNF6 (Santa Cruz Biotechnology, sc-13050X) 
antibody. A detailed protocol is provided as Supplemental Material. 

Generation of ChlP-seq libraries 

Sequencing libraries were generated from —10 ng of immunoprecipitated 
DNA using the Illumina TruSeq DNA sample preparation protocol with 
some changes: (1) Illumina TruSeq-indexed adapters were diluted to a 3:1 
adapter: insert ratio in each ligation reaction, (2) ligation products were 
size-selected (inserts of —130-280 bp) prior to PCR amplification, and (3) 
libraries generated from MNase-digested chromatin were size-selected 
after the amplification (inserts of a mononucleosome). 

Computational analysis 

Analysis details are available as Supplemental Material. Sequences were 
mapped to the mouse genome (version mm9) using bowtie (Langmead 
et al. 2009). Only those that mapped uniquely to the mouse genome were 
used for further analysis. A summary of the alignment results is provided 
in Supplemental Table 2. Analysis of rhythmic expression was performed 
as previously described (Menet et al. 2012). 

MNase-treated chromatin data sets Sequences were expanded to 
147 nucleotides (nt). Nucleosome signal was then retrieved at genomic 
locations of interest and normalized to the sequencing depth. These 
locations include CLOCK:BMALl binding (Koike et al. 2012), REV-ERBa 
binding (Feng et al. 2011; Cho et al. 2012), and DNase I-hypersensitive 
sites (Ling et al. 2010). The signal obtained for each library was averaged 
by time points. 

Sonicated chromatin data sets The CLOCK ChlP-seq signal was 
retrieved at genomic locations of interest, normalized to sequencing 
depth, and averaged for each library. 

Coassociation analysis The list of the publicly available mouse liver 
ChlP-seq data sets used in the coassociation analysis is provided as 
Supplemental Table 3 (see also the Supplemental Material). Coassociation 
analysis of transcription factor DNA-binding sites was performed using 
the algorithm described in the original ENCODE project paper (Dunham 
et al. 2012). The percentage of overlap (Fig. 4B), defined as the base-pair 
overlap ratio between two data sets, was calculated using the BEDTools 
suite (Quinlan and Hall 2010). 



Additional information 

All data sets are publicly available on Gene Expression Omnibus data- 
base at http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=vreftaseme- 
sumjw&acc=GSE4 7 1 45 . 
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